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Equivalence of Wilson Loops in AB JM and J\f = 4 SYM Theory 
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In previous investigations, it was found that four-sided polygonal light-like Wilson loops in ABJM 
theory calculated to two-loop order have the same form as the corresponding Wilson loop in TV = 4 
SYM at one-loop order. Here we study light-like polygonal Wilson loops with n cusps in planar 
three-dimensional Chern-Simons and ABJM theory to two loops. Remarkably, the result in ABJM 
theory precisely agrees with the corresponding Wilson loop in TV — 4 SYM at one-loop order for 
arbitrary n. In particular, anomalous conformal Ward identites allow for a so-called remainder 
function of conformal cross ratios for n > 6, which is found to be trivial at two loops in ABJM 
theory in the same way as it is trivial in TV = 4 SYM at one-loop order. Furthermore, the result for 
arbitrary n obtained here, allows for a further investigation of a Wilson loop / amplitude duality 
in ABJM theory, for which non-trivial evidence was recently found by a calculation of four-point 
amplitudes that match the Wilson loop in ABJM theory. 

PACS numbers: 11.25.Tq,11.15.Bt,11.15.Yc,11.30.Pb 



I. INTRODUCTION 

Our motivation to consider polygonal light-like Wilson 
loops in 3d Chern-Simons and TV — 6 superconformal 
Chern-Simons (ABJM) theory [l[ stems from the Wilson 
loop/scattering amplitude duality in TV = 4 super Yang- 
Mills. 

In planar TV = 4 super Yang-Mills n-particle MHV 
scattering amplitudes A^ HV = A^ ee M n are related to 
the expectation value of the n-cusped Wilson loop oper- 
ator 



(W n ) := -(0|TrPexp 



|0>. (1) 



The contour of the n-sided polygon C n is given by n 
points Xi (i — l,...,n) which are related to the massless 
particle momenta via Xi+i — x% = p». The segments of 



the contour are thus light-like, i.e. (Xi —Xi+i) 



Pi 



0. 



The relation between the Wilson loop and the scatter- 
ing amplitude is given by 



In M n = ln{W n ) + const . 



(2) 



This duality was discovered in the dual AdS$ x S 5 
string picture at strong gauge coupling in J2jand shown 
to exist also in the weak coupling regime [3j-l5| with pro- 
found consequences on the symmetries of these corre- 
lators leading to a dual superconformal [|| respectively 
Yangian symmetry 0] of scattering amplitudes, for re- 
views see [1,0. For a rieview of AdS/CFT integrability 
see [13 ■ 

It was found that the expectation value of the Wil- 
son loop in TV = 4 super Yang-Mills is governed by an 
anomalous conformal Ward identity that completely fixes 
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its form at 4 and 5 points and allows for an arbitrary func- 
tion of conformal invariants starting from 6 points. This 
so-called remainder function R n is indeed present start- 
ing from 6 points and leads to a correction [ll[ , [12] of the 
BDS ansatz [l3[ for planar gluon scattering amplitudes. 

Recently, the duality has been extended to amplitudes 
with arbitrary helicity states by introducing a suitable 
supersymmetric Wilson loop 0, . 

Furthermore, a duality between light-like Wilson loops 
with n cusps and n-point correlation functions of half- 
BPS protected operators in the limit where the positions 
of adjacent operators become light-like separated was es- 
tablished in [li-[20|. 

From the string perspective the scattering ampli- 
tude/Wilson loop duality in the AdS^/CFT^ system 
arises from a combination of bosonic and fermionic T- 
dualities under which the free AdS$ x S 5 superstring is 
self-dual [Hl[12|. Hence, for the existence of an analogue 
duality in ABJM theory one would require a similar self- 
duality of the AdS^ x CP3 superstring under the com- 
bined T-dualities. The problem was analysed in [23M28| 
but no T-self-duality could be established so far. 

At tree-level, recent developments have uncovered Yan- 
gian and dual superconformal symmetry of the ampli- 
tudes in ABJM theory [29I433]. In a vanishing result 
for the four-point one-loop amplitudes in ABJM theory 
was found and the authors speculated whether the two- 
loop scattering amplitudes in TV = 6 Chern-Simons could 
be simply related to the one-loop TV = 4 Yang- Mills am- 
plitudes. 

In (3f| we calculated the expectation value of the Wil- 
son loop operator JT]) in the planar limit for light-like 
polygonal contours C n in pure Chern-Simons and ABJM 
theory. 

Conformal Ward identities force (W n )i-i oop to depend 
only on conformally invariant cross ratios of the [xi — 
Xj) 2 . At one-loop order in pure Chern-Simons and ABJM 
theory we found that the correlators with four and six 
cusps vanish, leading to the conclusion that the allowed 
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function of conformal cross ratios is trivial at six points. 
Wc thus conjectured the n-point correlator to vanish 



loop 







(3) 



which was indeed proven in [361 ] and also non-trivial evi- 
dence for a duality between Wilson loops and correlators 
in ABJM theory was found at one-loop level. 

Furthermore, we computed the tetragonal Wilson loop 
Wi at two-loop order in pure Chern-Simons and ABJM 
theory. Remarkably, the result in dimensional reduc- 
tion rcgularisation with d = 3 — 2e for the correlator 
in ABJM theory is of the same functional form as the 
one-loop result in Af = 4 super Yang-Mills theory. Most 
interestingly, it was recently found by two independent 
approaches, using generalized unitarity methods in [37] 
and by a direct superspace Feynman diagram calculation 
in [38|, that the two- loop result for four-point scattering 
amplitudes in ABJM theory agrees with the Wilson loop 



Mj 2 > = 



A 



(2) 



A tree 
A 4 



— /Uf \ ABJM 



const. 



(4) 



upon a specific identification of the regularisation scales 
[37L [38j . This establishes the first non-trivial example 
for a Wilson loop / amplitude duality of the form in 
ABJM theory. Very recently, these result were extended 
to the more general case of ABJ theory in [39| . 

In light of these recent findings on structural similari- 
ties between observables in Af = 4 super Yang-Mills and 
ABJM theory, it is natural to ask, whether the duality 
between Wilson loops and amplitudes in ABJM theory 
continues to hold beyond n = 4, as it does in Af = 4. 

In this work we perform numerical computations to 
extend our findings of [H[ to the n-sided Wilson loop at 
two-loop order. Remarkably, we find that the hexagonal 
Wilson loop at two loops agrees with the corresponding 
Wilson loop in Af = 4 super Yang-Mills. 

We perform a detailed numerical analysis for the 
hexagonal Wilson loop leading to a guess for the n-point 
case, which we numerically check also for n > 6 in a lim- 
ited set of kinematical points. Again we find, that the 
result agrees with the result for the Wilson loop in Af = 4 
super Yang-Mills. It is thus natural to expect the result 
to hold for all n 



(wv ABJM 



/ 2- loop 
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where fi' — u 87re 7E , r„ is a constant that depends lin- 
early on n and is specified below (|19|) and where the finite 
contribution is given by the finite part of the Wil- 

son loop in Af = 4 SYM, which up to a constant 1 is the 



1 J"WL= J-BDS_„ /4((2) . 



finite part in the BDS conjecture i.e. for 
11 = 6 
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Thus, the Wilson loop in ABJM theory at two-loop 
order precisely agrees with the form of the Wilson loop 
or, via the amplitude / Wilson loop duality, with the 
n-point MHV amplitudes at one loop in Af = 4 super 
Yang-Mills. 

It would be very interesting to establish a six-point 
amplitude calculation in order to see whether the duality 
relation (J2J) in ABJM theory also holds at six points. Fur- 
thermore, it would be interesting to perform a four-point 
amplitude or Wilson loop computation at four loops, to 
check whether the conjectured BDS-like ansatz [38| for 
the four-point amplitude in ABJM indeed holds. 

The relations between Wilson loops, amplitudes and 
correlators seem to hold not only in the special case of 
the maximally supersymmetric Af = 4 SYM theory but 
also in ABJM theory. Since the duality may thus not 
just be a particular feature of Af = 4 SYM, it remains 
an important task to understand the precise origin of the 
similarity of these different observables in quantum field 
theory. 



II. N-SIDED WILSON LOOPS IN CS THEORY 

The solution of the Ward identity [35] for the light- 
like polygonal Wilson loop in pure Chern-Simons theory 
reads 



<W»>£Lp 



21n(2)]T ( 



2^2e 



2e 

9n{Uabcd) + 0(e) 



(7) 



where /2 2 = /x 2 7re 7E and g n is a function that de- 
pends only on conformally invariant cross ratios 

^abed — 

( x lb x ld)/( x ld x lb) which can be constructed 
starting from n = 6. 

At two loops there are two types of contributions 
to the Wilson loop in pure Chern-Simons theory, one 
from two-gluon diagrams (W„) two " gluon = J2 hjkl an d an- 
other one from diagrams involving a three-gluon vertex 

(w n y eitex = ^2i ijk 



(W n ) cs 



2-loop 



(W n ) 



two-gluon 



(W n ) 



(8) 



The indices i, j, k, I denote the edges that the propagators 
attach to, see figured] and their expressions are given in 



k 

(a) 



j 

(b) 



FIG. 1: 1(a) Two-gluon diagrams Iijkl given in (|A5[) . 1(b) 
vertex diagrams Jy*, given in (|A7I) . Dashed lines represent 
one or more light-like distances. 



appendix I A 11 lA~2l Contributions from gauge- and ghost- 
loops cancel in dimensional reduction regularisation [40j , 
for more details see (35[. 

As explained in [35| , the vertex diagrams are divergent 
in the region of integration where all three propagators 
approach the same edge (all diagrams with more than 
one propagator on the same edge vanish identically due 
to the antisymmetry of the Levi-Civita symbol) , and we 
split them up as in (|A8[) 



(w n y 



{W n ) dlv + (W n ) 



finite 



(9) 



The divergent part can be calculated analytically 



FIG. 2: Numerical results for the hexagon n = 6: (a) = 

^two-glucm = ( Wn )fi»te and their gum (c) = ffn ( Uabcd ) i 

see pO|l . for the kinematical points given in IA 31 We omit the 
common factor j (x) • 



tion yields the same constant 2 

geiuabed) =C6-121n(2), c 6 = 5.57 ± 0.05 . (11) 

In figure [5] we show the results for the two-gluon and 
vertex contributions for a continuously deformed kine- 
matical configuration, generated as explained in app. 
IA 31 in order to illustrate, how the different contributions 
vary while their sum remains constant. 




Thus, the function g n in ([7]) is given by 

9n(Uabcd) = (^„} tW °- glUOn + (^„> finit0 • (10) 

We evaluate these contributions using a Mathematica 
program that generates all n-point diagrams, performs 
the index-contractions and numerically integrates the 
diagrams for randomly generated kinematical configura- 
tions, for more details see appendix [Al 

Hexagonal Wilson loop 

For the hexagonal Wilson loop we evaluated the 
contributions in (|10p for a large set of conformally 
equivalent and conformally non-equivalent kinematical 
configurations. 

Conformally equivalent configurations must yield the 
same result, since, by the anomalous conformal Ward 
identity, the expectation value is constrained to the form 
([7]), and thus the function g n depends only on confor- 
mally invariant quantities. 

It turns out, that even for kinematical configurations 
which are not conformally equivalent, the unknown func- 



Generalization to n Cusps 

It turns out that also for n > 6 the function of 
conformal cross ratios is just a constant, i.e. 

9n{u a bcd) = c n - 2n ln(2) . (12) 

In figure[3]we show the dependence of the numerical con- 
stant c„ on the number of cusps n up to n — 14. Clearly, 
the constant depends linearly on the number of cusps n. 
It seems reasonable to assume that this dependence holds 
for all n, i.e. 

c n = a + b ■ n , (13) 

which is the line shown in fig. [3] with the parameters 3 
a = 6.6 ± 0.1, b = -2.028 ± 0.025. Thus, we expect the 
Chern-Simons contribution to the n-point Wilson loop to 



2 The analytical term with ln(2) in 1111 , 1121 arises from the multi- 
plication of the analytically known divergent term with an 0(e) 
expansion of the prefactor, see 1A6I , 

3 We determine the constants in 1131 1 from the results at n = 4 and 
n = 6, since here we have the smallest number of integrals and 
thus the best numerical result. At n-points we have to evaluate 
2(3) +2(") two-gluon ltA4t and vertex integrals (lAlOt . 
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4 V k 
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2e 



2M2)E 



(14) 



where g n is the constant given by (I12p and (|13p . 



(a) 




(c) 



FIG. 3: This figure shows the dependence of the constant c n 
(|13p on the number of cusps n of the Wilson loop. 



FIG. 4: The diagrams involving the one-loop corrected gluon 
propagator have the same structure as the 1-loop diagrams in 
N = 4 SYM. Dashed lines labelled with P and Q represent 
one or more light-like distances. 



loops, the power of 1/x 2 is 1 — 2e here, as opposed to 
1 — e in the one-loop computation in TV = 4 SYM. Thus, 
it is clear that the results will be very similar to the 
expectation value of the Wilson loop in TV = 4 SYM. 

As in TV = 4 SYM we have three classes of diagrams 
shown in figureHJ Diagram 4(a) vanishes due to the light- 
likeness of the edges, whereas 4(b) yields a divergent, and 



4(c) yields a finite contribution. 



We have \3l 



III. ABJM THEORY 

In ABJM theory we have two gauge fields A^, A^ and 
use the Wilson loop operator proposed in [4l|, which is 
a linear combination of two Wilson loops, each with one 
of the gauge fields, see also (35|. The one- loop contri- 
butions are identical to pure Chern Simons theory and 
cancel each other 4 . Both two-loop contributions yield 
the same result and thus, it is sufficient to calculate the 
Wilson loop with the gauge field A^. In addition to the 
Chern-Simons contributions (|14[) we have contributions 
of bosonic and fermionic matter fields which appear in 
the one-loop corrected gluon-propagator 



<^>2-fo J op =(Wn)£ oop + (W n ) 



2-loop 



(15) 



The matter contribution is similar to the one in j\f = 4 
SYM Q , since the one-loop corrected propagator calcu- 
lated in d dimensions 133 is 5 



1 

'N 



2-d T 



2\d-2 



(16) 



which up to two small differences is the tree level Af = 4 
SYM gluon propagator. The first difference is a trivial 
prefactor, and the second is that since we are at two 



toes? = ( j ) ( + V 2 ) £ h 



where 



(17) 



,y> 2 I ft 2 yi 2 ^v>2 



( XijSiSj x i+l 



S 7 S i X - ■ i -| S n S <i X 



,j + l S i S j) 



d-2 



There are n divergent diagrams It+i of the type shown 
in fig. |4(b)l 



1 ( X 1.i+2) 



2c 



and the finite diagrams Itj with \i — j\ > 1, sec fig. 4(c' 
were solved in [H 



Li 2 (l 



Li 2 (1 - ax 2 i+l j+l ) 



Li 2 (l - axl J+1 ) + Li 2 (l - ax 2 +1J ) 



where 



X i jX l 



x i,j+l x i+l,. 



4 They are also zero seperately, as explained above. 

5 We drop the derivative term, it would not contribute to the 
gauge-invariant Wilson loop as explained in [35t , 



The sum over all finite diagrams Y17>j+i A; = -^n^ i s 
related to the well-known finite part of the BDS conjec- 
ture via J"WL = j-bds _ n£( 2 ), explicitly spelled out 
in JSt for n = 4 and n — 6. 
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Then, the full matter part reads 



1. Two Gluon Diagrams 



(W n ) 



matter 
2-loop 



1 / N 



E 

i=l 



2e 2 



•4F. 



■WL 



(18) 



where we have restored the regularisation scale /i. Taking 
into account the Chern-Simons result (|T4|) . the full result 
in ABJM theory can be written as 



(W \ ABJM 



1 " (-x^ +2 M 2 87re^)^ 



2 

+ T 



■WL 



{2ef 



n 



(19) 



where r„ = -(f vr 2 +c„ - 2raln(2) - 5nln 2 (2))/4 and c„ 
is the numerical constant given by (|13[) . 

Indeed, this is of the same form as the one-loop result 
for the Wilson loop in M = 4 SYM. 



The contributions from the two-gluon diagrams are all 
finite and can be written as, see [351 ] . 

(W n }^°-^°»=Uj) £ (Ii^i+Iujk) (A4) 

^ ' i>j>k>l 

where 

J ({Zi - Zj) 2 ) {{z k - zi) 2 ) 

and the integration boundaries haven to be chosen ac- 
cording to the path ordering, such that z(sj) > z(sj) > 
z(s k ) > z(s{). The diagram vanishes due to the epsilon 
tensor contractions if the gluon propagator of at least one 
of the index pairs connects the same or adjacent edges. 

As a check of the numerics one can use the factorizing 
diagrams J"i+5,-i+3,i+2,i , which are just a p roduct of the 
analytically known one-loop diagrams [36l ] . 
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Appendix A: Two Loop Diagrams 

We use the same conventions as in [35J, i.e. the met- 
ric with rj^ v — diag(l, — 1, —1) and we define an n-sided 
polygon by n points Xi (i = l,...,n), with the edge i 
being the line connecting Xi and Xi+\. Defining 



n u n 

Pi = X i+1 - X i 



(Al) 

and parametrising the position z^ on edge i with the 
parameter Si 6 [0, 1] we have 

*g (A2) 
Furthermore, we use the notations 

e(p, q, r) = e fJ , up p' J/ q u r p and s t = 1 - s % . (A3) 
We use the Lagrangian of ABJM theory given in [35[ . 



2. Vertex Diagrams 

The contribution from the vertex diagrams is calcu- 
lated in the same way as in (35j . but generalized to n 
points and we find 6 



1 / N 



4 V k 



(7re 7E ) 2e (l-2e) (A6) 



I E ($+$) + °< e 



i>j>k 

The indices k indicate the edges the gluon- 
propagators connect to and 



p 



4jk = J ds i ^ }k ,d[l3} 3 e a0 ' < e tiaa e l/ i3xep' f rPtp' J j 

[d$Ad°dZA + d°Ad$dlA + dlAd?d$A] A 1 -' 1 



[dfAdfAdlA] A- d (l - d) 
where d!f = d/dzi^ and 



(A7) 



dffla= / dPi J , k 6(J20 m -l){p i j k )'V 2 - 2 , 

A = - zj k PA - z%Mi ■ 

Integrals with at least two propagators on the same edge 
(two identical indices i,j,k) vanish due to the antisym- 
metry of Iij k in the indices. Integrals I^l of the type 



Details will be presented in [42 
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i = k + 2, j = k + 1, k are divergent, see [35|, all other 
integrals are finite. The divergent integrals can be split 
up into a divergent and a finite piece 



(W n ) vertex = (W n f lv + (W n ) 



\ finite 



(A8) 



where the divergent piece can be evaluated analytically 
and reads 



(W r . 



\div 



\ 2 ^ 

\ i=l 

The reamining finite piece is 
(W n f nito 



N\ 2 (\n{2)^{~xl l+2 ^e^f 
k 



2c 



V 7 \ i>j>k j 

where rO'/),/^) are given by the following expressions 
for the case i = k + 2, j = k + 1, k. All other cases are 
treated purely numerically starting from (|A7[) . i.e. in 



particular 



r(«) 

l ijk 



in these cases. 



For i = 3,j = 2,k — 1, after solving two integrations 
and changing integration parameters to /3; = xy, f3j — 
xy, Pk = y (where x = 1 — x, y = 1 — y), we have 



'321 



In e 



1 t A A A (^) 1/2 



' 24 



' 24 



+ hi e' + /' 



13 



'13 



II 



1 (xxv) 1/2 
d Sj dxdy { »> (1 + \n{xxy) + ln(e'f')) 
e 7 



where the last term does not depend on the kinemati- 
cal quantities and can be further integrated analytically 
and/or evaluated to high numerical precision. We used 



the abbreviations 

f = -xlif = - x li x (y + ysjx) 

9 = ( x 13 + x 2i — X 14) X V ■ 

Furthermore, after solving two integrations we find 

These are the expressions we use for the numerical eval- 
uation of the finite parts of the divergent diagrams. 



3. Generation of kinematical configurations 

A set of n light-like (pj — 0) vectors p% satisfying mo- 
mentum conservation JZpf = can easily be generated 
by choosing the components of n — 3 vectors and the 
angle 9i between p\ — p° cos Q%, p\ — p\ sin 6{ . The re- 
maining components are then fixed. For n even it is 
possible to choose configurations, where all non-light-likc 



distances xfj are space-like. For the numerical evaluation 
we make use of this type of configurations, such that all 
integrals are real. 

The configurations used for the results shown in fig. 
[5] are obtained by continuously deforming two angles 
leading to conformally non-equivalent kinematical con- 
figurations. We use the angles 



.16 13 5 13 19 
i(a) := 7r < — a, — , —a, — , — . 1 

w ' 9 9 3 8 14' 



(All) 



and choose pf _0 = {1, —3, 4}, the remaining components 
are then fixed. The parameter a is chosen between a = 1 
and a = 1.2 in steps of 0.01. 
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